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Abstract 

Control and state estimation of nonlinear systems satisfying a Lipschitz continuity condition have 
been important topics in nonlinear system theory for over three decades, resulting in a substantial amount 
of hterature. The main criticism behind this approach, however, has been the restrictive nature of the 
Lipschitz continuity condition and the conservativeness of the related results. This work deals with an 
extension to this problem by introducing a more general family of nonlinear functions, namely one- 
sided Lipschitz functions. The corresponding class of systems is a superset of its well-known Lipschitz 
counterpart and possesses inherent advantages with respect to conservativeness. In this paper, first the 
problem of state observer design for this class of systems is established, the challenges are discussed and 
some analysis-oriented tools are provided. Then, a solution to the observer design problem is proposed 
in terms of nonlinear matrix inequalities which in turn are converted into numerically efficiently solvable 
linear matrix inequalities. 

Keywords: one-sided Lipschitz systems, quadratic inner-boundedness, nonlinear observers, 
nonlinear matrix inequalities, linear matrix inequalities 

I. INTRODUCTION 

THE observer design problem for nonlinear systems satisfying a Lipschitz continuity con- 
dition has been a topic of a constant research for the last three decades. Observers for 
Lipschitz systems were first considered by Thau in his seminal paper [1] where he obtained 
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a sufficient condition to ensure asymptotic stability of the observer error. Thau's condition 
provides a useful analysis tool but does not address the fundamental design problem. Inspired by 
Thau's work, several authors have studied observer design for Lipschitz systems using various 
approaches [|2|-[11|. Lipschitz systems constitute an important class of nonlinear systems for 



which observer design can be carried out using pseudo-linear techniques. Indeed, the Lipschitz 
continuity assumption provides a norm-based form of a nonlinear inequality which can be 
substituted into the observer error dynamics and, with some elaboration, the observer error 
dynamics can be cast in a numerically tractable format that is dominated by a linear term. 

The study of nonlinear Lipschitz systems is, however, overshadowed by some controversy 
surrounding the possibly restrictive nature of the Lipschitz continuity condition. The statement 
is somewhat controversial and, in fact, different researchers often refer to this condition as "very 
general and satisfied by most dynamical systems," and as "too restrictive to be of any practical 
use." There is, of course, some truth in both of these opposing views: most nonlinear equations 
satisfy a Lipschitz continuity condition, with the same Lipschitz constant, only locally; possibly 
in a small region of the state space. However, the Lipschitz condition is almost never satisfied 
globally thus placing a limit to the application of any result based on a Lipschitz continuity 
condition. 

The Lipschitz constant of such functions is usually region-based and often dramatically in- 
creases as the operating region is enlarged. On the other hand, even if the nonlinear system 
is Lipschitz in the region of interest, it is generally the case that the available observer design 
techniques can only stabilize the error dynamics for dynamical systems with small Lipschitz 
constants but, as discussed later, fails to provide a solution when the Lipschitz constant becomes 
large. The problem becomes worse when dealing with stijf systems. Stiffness means that the 
ordinary differential equation (ODE) admits a smooth solution with moderate derivatives, together 



with nonsmooth ("transient") solutions rapidly converging towards the smooth ones [12, p. 71]. 
In particular, piece-wise functions such as saturation, deadband and backlash which are common 
in control systems (e.g. in actuators) are known to impose stiffness into the underlying system of 
ODEs. This problem has been recognized in the mathematical literature and specially in the field 
of numerical analysis for some time and a powerful tool has developed to overcome this problem. 
This tool is a generalization of the Lipschitz continuity to a less restrictive condition known as 
one-sided Lipschitz continuity which has become one of the building blocks in numerical analysis 
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and has been extensively applied to the stability analysis of ODE solvers [13|-[15|. 

Inspired by the these advances in the mathematical literature, in this paper, we extend this 
concept to the nonlinear observer design problem and consider stabilization of the observer 
error dynamics based on the one-sided Lipschitz condition. The advantages gained through this 
approach are two-fold: i) Generalization: we will show that the one-sided Lipschitz continuity 
covers a broad family on nonlinear systems which includes the well-known Lipschitz systems 
as a special case, ii) Reduced conservativism: Observer design techniques based on Lipschitz 
functions can guarantee stability only for small values of Lipschitz constants which directly 
translates into small stability regions. All available results on Lipschitz systems, however, provide 
only sufficient conditions for stability and the actual observer might still work with larger 
Lipschitz constants, even though the tool used in the analysis and design are unable to provide 
theoretical evidence. The implication is that there is a significant degree of conservativeness in 
the Lipschitz formulation, a critique that has often been reported by researchers, but that has 
been difficult to correct and has produced no valuable alternative. 

In this work we provide this valuable alternative in the form of the one-sided Lipschitz 
condition. We will show that the one-sided Lipschitz condition generalizes the classical Lipschitz 
theory in the following sense: any dynamical system satisfying a Lipschitz condition satisfies 
also a one-sided Lipschitz condition. However, the one-sided Lipschitz constant is always smaller 
than its Lipschitz counterpart, a difference that can be significant even for very simple nonlinear 
functions p3|-p3|. Examples are presented illustrating this property as well as showing cases 
where a dynamical system satisfies a one-sided Lipschitz condition even-though it is not Lipschitz 
in the classical sense. Specially, when a dynamical system is stiff, the conventional Lipschitz 
constant inevitably becomes very large while one-sided Lipschitz constant is still moderate p3|- 
p3| . As a result, more efficient and less conservative observers can be developed in this context. 

Unfortunately, these major advantages come along with a greater degree of difficulty en- 
countered dealing with one-sided Lipschitz systems. Unlike Lipschitz functions, which lead to 
an inequality in a rather simple quadratic form, the one-sided Lipschitz formulation leads to 
a weighted bilinear form which imposes significant challenges in manipulating the Lyapunov 
derivative. 

Very recently, nonlinear observers for one-sided Lipschitz systems were considered in fT6l|, 
p7| . These references consider a problem similar in spirit to Thau's original work on Lipschitz 
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systems: given a one-sided Lipschitz system with a known Lipschitz constant and a pre-defined 
Lyapunov function candidate, obtain sufficient conditions for the existence of an observer. As 



with Thau's work, [ 16 1, [ 17 1 provide an important analysis tool but the (difficult) design problem 
remains unsolved. 

In this paper our goal is to acknowledge the advantages of the one-sided Lipschitz formulation 
over the conventional Lipschitz assumption in the control and observation theory, and in particular 
to formulate the observer design problem based on that. In this respect, not only do we provide 
basic analysis tools but also we address the design problem and present a complete solution. 
We introduce a new concept for nonlinear systems, the quadratic inner-boundedness, which is 
discussed in Section IV. We show that the class of quadratic inner-bounded systems also extends 
the class of Lipschitz systems. Our observer design solution is based on the one-sided Lipschitz 
property and the quadratic inner-boundedness of the system. 

The remainder of the paper is organized as follows: Section II introduces the one-sided 
Lipschitz condition and study its basis properties. In Section III we consider the observer problem 
based on this property and addressed observer stability. Section IV, which contains the main 
results, addresses observer design in the form of nonlinear matrix inequalities (NMIs). In the 
cycle, in order to use the efficient readily available numerical solvers, we convert the proposed 
NMI problem into linear matrix inequalities (LMIs). Section V presents an illustrative example. 

II. Mathematical Preliminaries and Problem Statement 

Throughout the paper M represents the field of real numbers, M" the set of n-tuples of real 
numbers and IR"^^ the set of real matrices of order nhy p. <, > is the (often called "natural") 
inner product in the space W\ i.e. given x,y E M", then < x,y >= x^y, where x'^ is the 
transpose of the (column vector) x G M". ||.|| is the vector 2-norm (the Euclidian norm) in M" 



defined by ||x|| = y/< x, x >. 

Consider now the following continuous-time nonlinear dynamical system 

x{t) = Ax{t) + u) Ae M"^" (1) 
y(t) = Cx{t) C e (2) 

where x G W^,u G W^,y G W and ^{x,u) represents a nonlinear function that is continuous 
with respect to both x and u. The system ([T])-([2]) is said to be locally Lipschitz in a region V 
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including the origin with respect to x, uniformly in u, if there exist a constant / > satisfying: 

-$(X2,M*)|| ^ -X2II y x,{t), X2{t) e V, (3) 

where u* is any admissible control signal. The smallest constant / > satisfying ([3]) is known 
as the Lipschitz constant. The region T) is the operational region or our region of interest. If the 
condition ([3]) is valid everywhere in M", then the function is said to be globally Lipschitz. 
The following definition introduces one-sided Lipschitz functions. 

Definition 1. [13| The nonlinear function is said to be one-sided Lipschitz if there 

exist p G M such that 

{^{xi,u*) - <^{x2,u*),xi - X2) ^p||xi-X2||^ Va;i,X2eP, (4) 

where p G M is called the one-sided Lipschitz constant. As in the case of Lipschitz functions, 
the smallest p satisfying (|4]) is called the one-sided Lipschitz constant. 

Similarly to the Lipschitz property, the one-sided Lipschitz property might be local or global. 
Note that while the Lipschitz constant must be positive, the one-sided Lipschitz constant can be 
positive, zero or even negative. For any function we have: 

I {^{xi,U*) - ^{X2,U*),XI- X2) I ^ \\^{xi,U*) - <I>(X2,U*)||||X1 -X2II 

and if ^{x,u) is Lipschitz, then: ^ l\\xi — X2\\^ ■ (5) 

Therefore, any Lipschitz function is also one-sided Lipschitz. The converse, however, is not true. 
For Lipschitz functions, 

- X2||^ ^ - $(x2,'U*),xi - X2) ^ /||a;i - X2||^, (6) 

which is a two-sided inequality v.s. the one-sided inequality in (|4]). If the nonlinear function 
m) satisfies the one-sided Lipschitz continuity condition globally in M", then the results 
are valid globally. For continuously differentiable nonlinear functions it is well-known that the 
smallest possible constant satisfying ([3]) {i.e., the Lipschitz constant) is the supremum of the 
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norm of Jacobian of the function over the region V (see for example p8|), that is: 

9$ 



/ = lim sup 



dx 



(7) 



Alternatively, the one-sided Lipschitz constant is associated with the logarithmic matrix norm 
(matrix measure) of the Jacobian. The logarithmic matrix norm of a matrix A is defined as [15]: 

\I + eA\\\ - 1 



lJi{A) = lim 



(8) 



where the symbol 



represents any matrix norm. Then, we have [15| 



p = lim sup 



dx 



Wx e V. 



(9) 



If the norm used in ([8]) is indeed the induced 2-norm (the spectral norm) then it can be shown that 

fi{A) = A, 



[19|. On the other hand, from the Fan's theorem (see for example [20|) 
we know that for any matrix, Xmax — '^rnax{A) = \\A\\ [20|. Therefore p < I. Usually 

one-sided Lipschitz constant can be found to be much smaller than the Lipschitz constant [ jl5| . 



Moreover, it is well-known in numerical analysis that for stiff ODE systems, p « I [14|, [15|. 
The one-sided Lipschitz continuity plays a vital role in solving such initial value problems. For 
instance, consider the following example which is taken from [[14} p. 173]. 



Example 1. Suppose ^(x) = —x^ with x G M. For the Lipschitz constant we have 

— $(X2)|| = \xl — xl\ < \xl + X1X2 + X2|.|Xi — X2I. 

By considering the line Xi = X2, we deduce that the Lipschitz constant on any set 

V = {x eR:\x\< r}, 

is necessarily ^ 3r^. Therefore, this function is not globally Lipschitz but is locally Lipschitz 
with a region-based Lipschitz constant which dramatically increases with the region of interest. 
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On the other hand, for the one-sided Lipschitz constant we have 

($(a;i) - $(x2), xi - X2) = -{x\- xl){xi - X2) 

= -{Xi - X2)'^{xl + X1X2 + xl) < -^(Xi - X2 f{xl + x1) < 0, 

which means that the function is globally one-sided Lipschitz with one-sided Lipschitz constant 
zero. We shall now illustrate an example which satisfies one-sided Lipschitz continuity but not 
Lipschitz continuity. 

Example 2. A simple example of a one-sided Lipschitz function which indeed is not Lipschitz 
is $(a;) = —sgn{x)^/\x\ with s G M, where sgn{.) denotes the sign (signum) function. This 
function is monotone decreasing and so globally one-sided Lipschitz with one-sided Lipschitz 
constant p = 0. In any interval x E [—m,m], the one-sided Lipschitz constant is [21]. 
On the other hand, this function is not Lipschitz in any domain containing the origin. Indeed, 
for X > 0, 

l/M-/(0)L^, 1 (10, 

\x-0\ \x\ ^\x\ 

which is unbounded as x — )■ 0. 



It is worth mentioning that the Lipschitz continuity property lies between continuity and 
continuous differentiability i.e. every Lipschitz function is continuous but not necessarily con- 
tinuously differentiable, while a one-sided Lipschitz function may be discontinuous. However, 
as in the case of Lipschitz systems, one-sided Lipschitz systems have only one solution (if 
exists) associated with each initial condition. Therefore, if the one-sided Lipschitz function is 
also continuous (a sufficient condition for the existence of a solution in ODEs), the system has a 



unique solution for every initial condition [22 p. 139]. It can also be shown that if a dynamical 
system is locally Lipschitz and globally one-sided Lipschitz with a strictly negative one-sided 



Lipschitz constant, then the system has a unique exponentially attractive equilibrium point [14|. 
A comprehensive treatment of the interesting properties of the one-sided Lipschitz functions is 
beyond the scope of this paper. The interested reader can consult the literature on numerical 



analysis, such as |14|, [13| and [15| and the references therein for further details. 
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III. Observer Structure 

In this section we consider the observation problem; i.e. given the dynamical system ([T])-([2]) we 
assume that only the input u and output y are available and study the feasibility of reconstructing 
the state x. To this end we consider an observer of the following form 

£{t) = Ax{t) + ^{x,u)+L{y-Cx). (11) 

The observer error dynamics is given by 

e{t) = x{t) - x{t), (12) 
e{t) = {A - LC)e{t) + ^ - ^ (13) 

where $ = ^{x,u) and $ = n). Using a Lyapunov function candidate of the form V(t) = 
e^Pe, the derivatives along the error trajectories are: 

V = e^[(A - LCfP + P{A - LC)]e + 2e^P($ - <l). (14) 

Our goal is, assuming that $ satisfies a one-sided Lipschitz condition, find an observer gain 



L such that the observer error dynamics is asymptotically stable. Equation (14), however, is 
independent of the assumption made on the function $ and can be used when $ is either 
Lipschitz or one-sided Lipschitz. Before further exploring the one-sided Lipschitz problem we 
digress to compare with the Lipschitz case. 

A. Motivating problem 



Given the form of equation (14), the standard approach used in the Lipschitz case is to find L 
such that {A - LCY P + P{A - LC) = -Q is negative definite and e^Qe > 2e^P($-<l>). This 
approach, first proposed by Thau in 1973 [[TJ, has dominated the literature on Lipschitz systems 
ever since. The implicit idea behind this approach is to use of the output injection term in the 
observer dynamics to ensure that the linear part of the observer error dominates the nonlinear 
terms. This, in turn, is facilitated by the strong square norm condition ([3]) satisfied by Lipschitz 
systems which leads to the conservative nature of the result. In other words, in the process of 



employing the Lipschitz property (|3]), the term e^P($ — $) in the Lyapunov derivative ([14]) is 
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replaced by a strictly positive term forcing the rest of the derivative to be sufficiently negative 
to compensate the remaining terms. 

It is important to note that the term e^[(A - LCY P + P{A - LC)]e < if and only 
if y4 — LC has eigenvalues with negative real part. Unlike the Lipschitz constant I, which is 
positive by definition, the one-sided Lipschitz constant p can be any real number. Thus, the term 
2e^P($ — $) in ( [141 ) '^^^ be negative. Hence, a negative Lyapunov derivative may be guaranteed 
even with a positive definite {A — LC)^ P + P{A — LC) and consequently the linear terms in 
A — LC is not necessarily required to have eigenvalues with negative real part. This means that 
the linear terms not necessarily dominate the nonlinear function $, which in turn can lead to 
less conservative results. 

The mathematical description behind the need of A — LC having eigenvalues in the left half 
plane can be traced back through a substantial body of literature for Lipschitz systems such as 
[[T|, [|2|, ||5|-|(8| while the freedom of the one-sided counterpart from such necessity is established 
in this article. 

B. Challenging Obstacle 



We now return to our main objective and endeavor to find L that makes F < in ( 14). This 
problem is nontrivial. Unfortunately, to the best of the authors' knowledge, there is no result in 
mathematics relating the weighted bilinear form 2e^P($ — $) to the bilinear form 2e^($ — $) 
which is required to take full advantage of the properties offered by the one-sided Lipschitz 
formulation. This is the main obstacle to overcome in this work. A very simple first approach 



to this problem is to consider P = L With this choice, from (14) we have: 

V = e^[{A - LCf + {A- LC)]e + 2e^(<l> - <l>) < e^[{A - LCf + {A- LC) + 2p/]e. 
where we substituted e^{^ — < pe^e. Hence, in order to have y < 0, we must have 

{A - LC)'^ + {A- LC) + 2p/ < ^ p{A - LC) < -p. (15) 



Inequality ( [15] ) is an LMI which can be efficiently solved using any available LMI solver to find 



the observer gain L. For the logarithmic matrix norm the following inequality can be used [19]: 

-^l{-A) <mi{A) < ii{A), z = l,...,n. (16) 
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Therefore, - LC)) < ^\i{A - LC). On the other hand, we want ii{A - LC) < -p, 

so as a necessary condition, we must have 

max3fJAi(A - LC) < -p. (17) 

i 

Furthermore, suppose A — LC is not stable (not stabilizable). We can always find a > such that 
(A — LC — al) is stable where a > maxj ^Xi{A — LC). Then the observer error is asymptotically 
stable if 

{A - LCf + {A- LC) + 2pl ={A-LC- al)'^ + {A-LC- al) + 2al + 2pJ < (18) 



Now, a sufficient condition for ( 1 8 1 to be true is 



=^a + p<Q^p<-a< -max3?Ai(A - LC). 

i 

Remark 1. In general, defining P = I h too restrictive to be regarded as a viable solution as 
it forces the error dynamics to be stable with the Lyapunov candidate V = e^e. Alternatively, 
we can redefine the one-sided Lipschitz property as 

{P^{Xi,U*)- P^{X2,U*),XI-X2) ^ p\\x,-X2f V Xi (t) , X2 (t) G P. (19) 

With this definition, we obtain e^P($ — $) < pe^e and the resulting stability condition is 

{A - LC)'^P + P{A - LC) + 2pl < 0, (20) 



which is an LMI. The idea is to find p, P and L such that ( [T9| ) and ( [20| ) are simultaneously 



satisfied. While ( [20| ) is a definite programming, ( [191 ) is a function inequality that needs to be 
satisfied for all points inside the region V for all times. 

Solving the two inequalities together is also an open problem, much more difficult than finding 
p in the standard one-sided definition (|4]) which is still a challenging task. An alternative is to first 



find P and p from ([19]) (which is also very difficult) and then solve (20) for L. This approach. 



however, will not bring additional benefits as it is as restrictive as the previous solution in the 



sense that it forces the error dynamics to be stable with the specific P known apriori from (19) 



This is similar to the analysis result of [[T6|, [17| in which P is assumed to be known apriori. 



The above discussion provides some analysis insight but does not address the fundamental 
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design problem in a satisfactory manner. In the next section we propose a complete solution to 
this rather involved design problem. 

IV. Main Results 

In this section, we first introduce the concept of quadratic inner-boundedness for the function 
m). Our design solution will make extensive use of this concept. 



Definition 2. The nonlinear function $(a;,M) is called quadratically inner-bounded in the 
region ^ if V xi, X2 G -D there exist /3, 7 G M such that 

- ^{X2,u)f{^{xi,u) - ^{X2,U)) < (3\\Xi -Xsf +7 (Xi - X2,^{Xi,u) - $(x2,m)) • 

(21) 

It is clear that any Lipschitz function is also quadratically inner-bounded (e.g. with 7 = 
and (3 > 0). Thus, Lipschitz continuity implies quadratic inner-boundedness. The converse is. 



however, not true. We emphasize that 7 in pT] ) can be any real number and is not necessarily 
positive. In fact, if 7 is restricted to be positive, then from the above definition, it can be concluded 
that $ must be Lipschitz which is only a special case of our proposed class of systems. We will 
discuss the special Lipschitz cases in details, later in this section. Figure [T] shows the relation 
between the Lipschitz, one-sided Lipschitz and quadratically inner-bounded function sets. 




Fig. 1. The Lipschitz, one-sided Lipschitz and quadratically inner-bounded function sets 



From now on we assume that is one-sided Lipschitz in V and quadratically inner- 

bounded in v. All of our results will be valid in the intersection VnV (the operational region). 
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If is both globally one-sided Lipschitz and globally quadratically inner-bounded, then 

all of our results will be valid globally. With the above notation, the following inequality holds 
for the estimation error. 



(22) 



In the following Theorem, we propose a method for observer design for one-sided Lipschitz 
systems. 

Theorem 1. Consider a nonlinear system satisfying inequalities Q and pT] ) with constants p. 



/3 and 7, along with the observer ( [TT] ). The observer error dynamics is asymptotically stable if 
there exist positive definite matrix P, symmetric matrix Q, matrix L and a positive scalar a > 
such that the following matrix inequalities problem is feasible: 



{A - LCfP + P{A - LC) < -g, 



7 + 2a > 
A 



1) < 



where ^ = {(3 + 1) + p{j + 2a). 



Proof: We have 



(23) 
(24) 
(25) 
(26) 



V = e^[{A - LCfP + P{A - LC)]e + 2e'^P($ - <l). 



(27) 



for any nonzero a we can write 



1 r " 1 r 

2e^P(<l>-<l>) = - e + «($-$) P e + «($-$) 

--e^Pe-a(<!>-^fP(<i>-^), Va G M - {0}. 
a 



(28) 



Assuming a > 0, then 

2e^P(<l> - !>) < -A„„^(P) [e + a(<l> - <l)l ^ [e + a(<l> - <I>)1 - -e^Pe - a(<l> - <|))'^P($ - 

(29) 
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Using the quadratic inner-boundedness property we have 



e^e + a^{^ - <l')^($ - $) < (/3 + l)e^e + 7e'^($ - $) + (a^ - 1)(<I> - <l)^($ - <l) 

^ [e + a(<l> - <l)j^ [e + a($ - <l)j < (/3 + l)e^e + (7 + 2a)e^($ - $) + (a^ - 1)||$ - <lf . 

(30) 



Substituting (30) into (29) leads to 



2e^P($ - $) <-A^,,(P) (/3 + l)e^e + (7 + 2a)e^($ - $) + {a^ - 1)||$ - $| 
a 



a 



-e^Pe-«(<l>-<l>)^P($-$). 



(31) 



Based on the Rayleigh's inequality, for any a; > we have 



($ _ p($ - $) > A™„(p)ii<i> - - -^)< -aXminipm - •^'ii- 



Hence, from ( [3T] ), using the one-sided Lipschitz inequality (|4]) and knowing that 7 + 2q; > 0, 
we obtain 



2e^P(<l> -<!>)< 



a 



1$ _ <|)||2 



+ -A™„.(P) + 1) + p(7 + 2«)] e^e - -e^Pe. 



(32) 



We know that K(P)(a^ ~ 1) < Q^^^ where k(P) is the condition number of P or ^Ama2.(P)(a^ — 
(P) < 0. Then, 

2e^P($ - $) < -A^,,(P) + 1) + p(7 + 2a)] e^e - -e^Pe. (33) 

a a 



Now we substitute (33) into the Lyapunov derivative. We obtain 



V = e' [{A - LCy P + P{A- LC)]e + 2e' P($ - $) 

< \(A - LCYP + PiA - LC)] e + -A„„^(P)e^e - -e^Pe, (34) 
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where ^ = (/3 + 1) + p(7 + 2a). Therefore, in order to have V < it suffices to have: 



{A - LCfP + P{A - LC) + ^\ma.{P)I - -P < 

a a 



a 



{A - LCfP + P{A - LC) - -P 



a 



< 0. 



(35) 



For any two symmetric matrices A and B, it can be shown that Aj < \i{A) + \i{B), where AjS 



are the sorted eigenvalues |20|. Thus, 



A. 



{A - LC)^P + P{A - LC) -^P < ^max [{A - LC)^P + P(A - LC)] + X^^x (^-0 



Xmax [{A - LC)^P + P{A - LC)] - -XrmniP). 

a 



(36) 



Now without loss of generality, we assume that there exists a symmetric matrix Q, such that 
{A — LC)^P + P{A — LC) < —Q. Note that Q is not necessarily positive definite (meaning 
that (A — LC) is not necessarily stable). Thus, 



X^ax [{A - LC)^P + P{A - LC)] < Xn,ax{-Q) = -XminiQ). 



(37) 



Substituting from (37), (36) and (35) into (34), we get 



^ 1 

Xmaxi^P) XfniniyP) Xmini^^) ^ 0. 

a a 



This completes the proof. 



(38) 



Remark 2. From a second glance at ( [26] ), it is clear that if — 1 < (0 < a < 1, taking 
a > into consideration), then the inequality is always satisfied. On the other hand, if — 1 > 
(a > 1, knowing that a > 0), we have k(P) < Then, as a approached one, the admissible 



condition number of P becomes larger. If 7 > —2, we can always pick a = 1 so that both (25) 



and ( [26] ) are always satisfied. On the other hand, if ^ < 0, then ( [24] ) is always satisfied. So, if 
in addition, a can be chosen such that < 0, then the NMIs in Theorem 1 reduce to a simple 



LMI feasibility problem (solving only (23)) 
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A. Special Lipschitz cases 

As mentioned earlier, the class of nonlinear systems we are addressing here covers Lipschitz 
systems as a special case. It is worth to briefly study under what conditions one-sided Lipschitz 
systems reduce to the Lipschitz case. 



• Case 1. If pT] ) is satisfied with 7 = 0, then (3 has to be positive and $(a;,t) is Lipschitz 
with Lipschitz constant ^/]3. 

• Case 2. If for a one-sided Lipschitz nonlinear function <l>(x, u), pT) is satisfied with 7 > 0, 
then l3 + 'yp has to be positive and it is easy to show that $(a;, u) is Lipschitz with Lipschitz 
constant a//3 + 7p. 

• Case 3. As mentioned in Remark 2, if — 1 < (0 < a < 1), then the inequality ( [261 ) 
always satisfied. On the other hand, from ( [SO] ) and knowing that — 1 < and 7 + 2a; > 0, 
we can write 

r " 1 ^ r -^1 

e + «($-$) e + a($-$) < (/3 + l)e^e + (7 + 2a)e'^($ - $) 
< + 1) + (7 + 2a)p]e^e \\e + «($ - $) f < [(/? + 1) + (7 + 2a)p] ||ef . (39) 



As a consequence, we have 



|e|| + a||<l> - <l|| < ||e + «($ - l>)|| < A/(/3Tiy+77T2ajp ||e| 



1$ - $11 < 



a 



1 + v/(/3 + 1) + (7 + 2a)p 



(40) 



which means that $ is Lipschitz. In other words, only for Lipschitz functions we can find 



< a < 1 such that 7 + 2a; > and (33) are both satisfied 



B. LMI formulation 

Theorem 1, provides a design method for nonlinear observers for one-sided Lipschitz systems 
in the form of the nonlinear matrix inequalities (NMIs) (|23])-([26l). The difficulty, however, is that 
although Theorem 1 provides a legitimate solution to our problem, there is currently no efficient 
solution in the numerical analysis literature capable of solving NMIs. Unlike the nonlinear case, 
however, linear matrix inequalities (LMIs) can be efficiently solved using commercially available 
packages such as the Matlab LMI solver. We now show how to cast the proposed nonlinear matrix 
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inequalities solution into the LMI framework to take advantage of the efficient numerical LMI 
solvers readily available. 

Using Fan's theorem mentioned earlier, we can write 

X^a.[{A - LCfP + P{A - LC)] < 2a^a4P{A - LC)] < 2X^a.{P)amax{A - LC). (41) 



Substituting this back to p5) yields to 

{A -LC) <0 

a a 



X'max i,P) 



2crmax{A - LC) + 



a 



< -XmUP) l^{P) 

a 



2arnax{A-LC) + 



< 



a 



(^max^A — LC) < 



(42) 



2a \k{P) 

which by means of Schur's complement and change of variable A = ^ is equivalent to the LMI 



( |43| ). LMI ( |45| ) represents the condition k(q;^ — 1) < 

Based on the above discussion which also serves as the proof, the following corollary provides 
an LMI solution to our observer design problem. 



Corollary 1. Consider a nonlinear system satisfying inequalities Q and pT) with constants 



p, (3 and 7, along with the observer (H) . The observer error dynamics is asymptotically stable 
if there exists a matrix L and positive scalars a > and < A < 1 such that the following 
matrix inequalities problem is feasible: 



^(A-0/ (A-LC)^ 
(A-LC) iiiX-Ol 

7 + 2a > 0, 
A > 1 - At, 



>0, 



where ^ = (/3 + 1) + p(7 + 2a). 



(43) 

(44) 
(45) 



We now summarize the observer design procedure based on Corollary 1. 
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C. Design Procedure 

• Step 1: Pick an a > such that 2a + 7 > 0. 

. Step 2: If p = and /3 > 0, Stop; otherwise, calculate ^ = /3 + p(7 + 2a) + 1. 

• Step 3: Check if the conditions A — ^ > and < A < 1 are consistent. If Yes, go to Step 
4; otherwise, adjust a (if p > then decrease a, otherwise increase a) and go to Step 1. 

• Step 4: Solve the LMIs in Corollary 1 for L and A. 

Note that if Step 3 is passed, the LMIs in Step 4 are always feasible meaning that an observer gain 
L will always be found. In order for a feasible solution to exist, we require both 2q! + 7 > and 
^ < 1 to be satisfied. Therefore, one can choose a > such that 2q;+7 > and ^+p(7+2q;) < 
are simultaneously satisfied. This guarantees the feasibility of the solution space in Corollary 
1. A structural limitation here is that if p = 0, then /3 has to be negative. The variable k>{— j) 
calculated in Corollary 1 is the condition number of the P matrix used in the Lyapunov function. 
Any P with such condition number would be acceptable. Although easy to do, as our goal of 
finding the observer gain L is already achieved, this step (finding P) is unnecessary except for 
analysis purposes of the results. 

Remark 3. The beauty of the above result is in its simplicity and in having the LMI form. 
It is, however, worth mentioning that this solution is conservative in the sense that we have 

replaced Xmax[{A — LC)^P + P{A — LC)] which may or may not be positive, with the positive 
quantity Xmax{P)o'maxiA — LC). Nonetheless, the solution provides a major improvement over 
the original solution in Theorem 1, in terms of numerical solvability. 

Remark 4. There is yet another advantage in the LMI formulation. The one-sided Lipschitz 
constant p appears linearly in the proposed LMIs. Hence, one can take advantage of the convexity 
of the solution space of the LMIs and solve an LMI optimization problem, maximizing the 
admissible one-sided Lipschitz constant for which the observer error stability is guaranteed. The 
resulting maximized p acts as an upper bound to those of nonlinear functions $(s). This is 
specially helpful in the situations that an accurate estimate of the one-sided Lipschitz constant 
is hard to obtain. Alternatively, the optimization can be performed over the quadratic inner- 
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boundedness parameters, ^ or 7. 

V. Illustrative Example 

In this section, we illustrate the proposed observer design procedure through a numerical 
example. 

Example 3. Suppose that the equations of motion of a moving object are given in the 2D 
Cartesian coordinates as follows: 

X — X — y — x{x^ + y^), 
y^x + y-y{x'^ + y'^), 



and y is measured. We define the state vector as x 



x y 



and y = y as the output, in 



which the variables are bolded to avoid ambiguity. We have: 



1 -1 
1 1 

1 



x + 



-x{x^ + y^) 

-y{x^ + y^) 



Then we verify the one-sided Lipschitz property. 
($(xi) - $(x2),xi - X2) = {X1X2 + ym) [{xl + yl) + {xl + yl)] - {xl + ylf - {xl + ylf 



On the other hand, XxX2 = \{xl-\- x\) — — ^2)^. Thus, 

X1X2 + ym = ^ [||xi||^ + ||x2||^ - ||xi - X2IP] . 



(46) 



Therefore, 



($(Xi) - $(X2),X1 - X2) = - [||Xi||2 + IIX2II' - ||Xi - X2II'] . [||Xi||2 + IIX2II'] - ||Xi| 



-^l|Xl-X2f [||Xif 



1 



< 0. 



X2 



(47) 



This means that the systems is globally one-sided Lipschitz with the one-sided Lipschitz constant 
p — Q. Now, lets verify the Lipschitz continuity property. $ is continuously differentiable, so an 
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estimate for the Lipschitz constant is the supremum of the norm of the Jacobian matrix, J. We 
have: 



J 



-3x^ — y"^ —2xy 
—2xy —3?/^ — 



J is a symmetric matrix; therefore, its induced 2-norm equals its spectral radius. 

\\J\\ = <yn.aAJ) = max |Ai(J)| = max(| - 2{x^ + y^) ± (x' + y^)\) = ^{x^ + y^) = 3r^ 

i 

This means that the system is locally Lipschitz and on any set 

I? = {x G : ||x|| < r}, 

the Lipschitz constant / is 3r^, i.e. the Lipschitz constant rapidly increases with the increase of 
r. We will come back to this later. Now we proceed our one-sided Lipschitz based observer 
design. We need to verify the quadratic inner-boundedness property of the system, as well. The 



left hand side of the quadratic inner-boundedness (21 ) is: 



LHS : [$(xi) - $(X2)]^ [$(xi) - $(X2)] 

= {x\ + ylf + {xl + ylf - 2(xiX2 + yiy2){x\ + yl){xl + yl) 
= \\xi\\^ + llxaf - [||xif + \\x2f - ||xi - Xaf ] ||xi||2||x2|p 

= llXif [ll^lf - M"] + WM" [M' - llXlf] + ||Xi - X2|n|Xif IIX2I 

ixiir - llx2lr]^ [ii^if + ii^^f] + iixi - X2f iix,inix2ir, 



in which (46) is used. The right hand side of the quadratic inner-boundedness (21) is: 



RHS : 7($(xi) - $ (X2), Xl - X2) + f3\\x^ - x^l 



-|||Xi-X2f [ll^if 



Xl - X2 



7 
2 



-||X2||2]' + /3||Xi-X2| 
Xlir-||X2f]', 



in which (47) is used. We have to find values for (3 and 7 and a region V such that for all 
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X eV, LHS < RHS. Comparing the two, it suffices to have: 

II ||2 I II ||2 7 

||Xi|| + IIX2II S — 



\x,r.\\x,r<^-i[\\x,r+\\x,r]'<^+^. 



Considering the set 



V = {xeR'^ : \\x\\ < r}, 



it suffices to have: 



- 2 - V 4 ' 



Hence, 




r = min| ^/-- , 1 , 7 < 0, /3 + ^ > 0. (48) 

Also, since p = 0, ^ = /3 + l and thus according the LMIs in Corollary 1, since A < 1, in order 
to have \ — ^ > 0, (3 has to be negative. As the system is globally one-sided Lipschitz (V = M?), 
V nV = v. It is clear that by choosing appropriate values for 7 and (3, the region V can be 
made arbitrarily large. If we take /3 = —200, = —199) and 7 = —141, we get r = 5.9372. 
Then we take a = 70.6 (to ensure 7 + 2a > 0) and solve the LMIs in Corollary 1. We get: 

1 T 



A = 0.999892, L 



-1.000000 1.000000 



(49) 



Figure (|2]) shows the system trajectories along with their estimates and the system phase plane. 
This system has a stable limit cycle at + = 1 which is inside V. As this limit cycle 
is a global attractor, the region V is the region of attraction for the observer error dynamics. 
For comparison purposes we now consider the conventional Lipschitz formulation. With r = 
5.9372, the corresponding Lipschitz constant is / = 3r^ = 105.75 (compare it with the one- 
sided Lipschitz constant p = 0). It is highly unlikely that an observer designed based on the 
conventional Lipschitz approach can work with such a large Lipschitz constant. The maximum 
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Fig. 2. The system and observer states and the phase plane 



Lipschitz constant that those observers can handle are normally at least an order of magnitude 
less than this. For example, the maximum admissible Lipschitz constant for the observer designed 
using [9] in this case is / = 1.0324. The advantage of the new approach in this case, is evident. 

1 

is indeed unstable, confirming our finding of A — LC not being 



1 



Note that A-LC = 

necessarily stable. To verify our result, lets calculate the Lyapunov derivative. Any positive 



definite matrix P with condition number k = \ h acceptable. Lets take P 



\ 
1 



. From 
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(|35l) we obtain 

^ -0.4187 < 0. (50) 



a 



{A - LCfP + P{A - LC) - -P 

a 



VI. Conclusions 

This article introduces a new class of nonlinear systems, based on the so-called one-sided 
Lipschitz continuity condition commonly used in numerical analysis, as a generalization of the 
well-known class of Lipschitz systems. The observer design problem for this class of systems, 
namely one-sided Lipschitz systems, is established. The advantages of designing observers in 
this context are explained and the challenges discussed. The mathematical preliminaries of the 
addressed class of systems are studied and some analysis-oriented tools for the observer stability 
are provided. Then the challenging design problem is tackled and a solution is proposed based 
on nonlinear matrix inequalities which are in the cycle casted into the linear matrix inequalities. 
In order to obtain the proposed solution, a new property for the nonlinear functions, the quadratic 
inner-boundedness is introduced. Finally, a straightforward observer design procedure is given 
that can be easily applied to the considered class of system using the available numerically 
efficient LMI solvers. The efficiency of the approach is shown through an illustrative example. 
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